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Parameterizing the Power Spectrum: 
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The power spectrum is traditionally parameterized by a truncated Taylor series: ln'P(fc) = In 7-* + 
(n* — 1) ln(fc/fc„) + |n', ln 2 (fc/fc„). It is reasonable to truncate the Taylor series if |n'„ ln(fe/fe*)| <C 
1 71* — lj, but it is not if |n» ln(fe/fc») | > |n* — lj. We argue that there is no good theoretical reason to 
prefer \n' t | <C \n* — 1|, and show that current observations are consistent with |n* ln(fe/fe*)| ~ |n* — 1| 
even for [fn(fe/fe*)| ~ 1. Thus, there are regions of parameter space, which are both theoretically 
and observationally relevant, for which the traditional truncated Taylor series parameterization 
is inconsistent, and hence it can lead to incorrect parameter estimations. Motivated by this, we 
propose a simple extension of the traditional parameterization, which uses no extra parameters, but 
that, unlike the traditional approach, covers well motivated inflationary spectra with |n'„| ~ |n* — 1|. 
Our parameterization therefore covers not only standard-slow-roll inflation models but also a much 
wider class of inflation models. We use this parameterization to perform a likelihood analysis for 
the cosmological parameters. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

In estimating the cosmological parameters and the pri- 
mordial power spectrum using the observational data, 
such as WMAP and SDSS [3, HQ, one introduces free 
parameters to parameterize the power spectrum V(k). 
There is some arbitrariness in how to do this beyond 
simply parameterizing the amplitude of a scale invari- 
ant spectrum. It has become traditional to consider a 
truncated Taylor expansion about some particular pivot 
scale fc*, including two derivative terms, the spectral in- 
dex n — 1 = dlnV/dhik and its running n' = dn/dhik. 
This traditional parameterization is motivated by sim- 
plicity and the standard slow-roll approximation. 

The standard slow-roll approximation is satisfied by 
the simplest single-component models of inflation, but 
one also has to be prepared for not so simple single- 
component models, and it may generically not be sat- 
isfied by multi-component models of inflation. Cur- 
rently there is no observational reason to disfavor multi- 
component models of inflation, and there are some the- 
oretical reasons to prefer them in terms of naturalness 
0. Thus using the standard slow-roll approximation to 
justify a parameterization of the power spectrum is dan- 
gerous. 

The truncated Taylor expansion could give a poor ap- 



proximation if the range of k of our interest extends far 
from fc*, unless we take a sufficient number of higher 
derivative terms in the Taylor expansion. However, in- 
creasing the number of derivative terms to as many as one 
desires would not be a practical approach for parameter 
estimation. 

In this paper, we propose an improved parameteriza- 
tion of the power spectrum, which is as simple as the 
traditional truncated Taylor expansion in that it uses 
the same number of parameters, reproduces the tradi- 
tional truncated Taylor expansion in the standard slow- 
roll limit, but has a better motivated form outside of this 
limit. Because our parameterization has the same num- 
ber of free parameters as the traditional truncated Tay- 
lor expansion, it can be straightforwardly implemented 
in existing numerical codes. 

The layout of our paper is as follows. The trun- 
cated Taylor expansion and its drawbacks are discussed 
in Sec. [H] Our proposed parameterization is presented 
in Sec. lIII Al and its inflationary motivation in Sec. lIII Bl 
In Sec. IIVI cosmological parameter estimation via the 
Markov-Chain Monte-Carlo (MCMC) method using our 
parameterization is compared with that using the tradi- 
tional truncated Taylor expansion. Conclusion and dis- 
cussion are in Sec. E] 



II. THE TRADITIONAL APPROACH 



*On sabbatical leave from Department of Physics, KAIST, Repub- 
lic of Korea. 



We describe the traditional approach, explain the 
weakness of its traditional justification, and describe 
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some possible undesirable consequences of adopting it. 

A. The traditional approach: a truncated Taylor 

expansion of the power spectrum 

Traditionally, the deviation from a scale-invariant 
spectrum is quantified by the tilt n — 1 = d In V/d In k 
and its running n' = dn/dhik by Taylor expanding in 
log-log space around a pivot point In fc* 

\nV(k) = lnP* + (n. - 1) In + ±< In 2 ^ 

"(2-1) 

This traditional truncated Taylor expansion approach as- 
sumes the second and higher derivatives of n are negli- 
gible, which is not a trivial assumption, especially for 
|m(fe/fe*)| > 1. Generally, one might expect this to 
be a good assumption if |ra* ln(fc/fc*)| <C |n* — 1|, but 
a bad one if ln(fc/fc*)| > \n* — 1|. We shall show 
in Sec. IIVI that current observations are consistent with 
ln(fc/fc*)| ~ |n* — 1|, and hence that there is no ob- 
servational justification for this assumption. 

B. The traditional justification of the traditional 
approach: the standard slow-roll approximation 



allows extra freedom to build more natural models, see 
for example Ref. [j]. In these multi-component models, 
the relevant slow-roll parameters 1 are small but, without 
fine-tuning, most of the irrelevant slow-roll parameters 
will not be small. The non-small irrelevant slow-roll pa- 
rameters then tend to cause the relevant slow-roll param- 
eters not to vary slowly. Thus, although standard slow- 
roll could be considered the generic, though not exclu- 
sive, case in single-component models, non-slowly vary- 
ing slow-roll parameters may be the more generic case 
in multi-component models. Thus, if one wants to cover 
all reasonable models and hence be able to distinguish 
amongst them, one should relax the assumption of slowly 
varying slow-roll parameters. 

The general slow-roll approximation 0, Q drops this 
extra assumption of slowly varying slow-roll parameters, 
covering the cases of 

|n-l|>|n'|>|n"|>... (2.4) 

Hence the general slow-roll approximation includes 
Eq. (|2.3[1 as a special case, so that we can test the as- 
sumption of the standard slow-roll approximation, rather 
than assuming it a priori. 



The standard slow-roll approximation assumes that 
the inflationary slow-roll parameters are both small and 
slowly varying. In terms of observable quantities, small- 
ness of the relevant parameters translates to 

|n-l|<l (2.2) 

which is required by observations, while the slowly vary- 
ing condition translates to the hierarchy 

(2.3) 

Thus we see that the validity of the traditional truncated 
Taylor series approach is equivalent to the assumption of 
slowly varying slow-roll parameters. Thus if we assume 
the standard slow-roll approximation, as is often done, 
neglection of the higher order derivatives in Eq. (|2.1|) 
follows. 

In terms of inflationary model building, for example 
to get enough inflation, there is no need for the slow-roll 
parameters to be either small or slowly varying. How- 
ever, as noted above, smallness of the relevant parame- 
ters is required for approximate scale invariance of the 
spectrum. In many simple single-component models of 
inflation, the requirement of small slow-roll parameters 
forces the slow-roll parameters to be slowly varying, but 
there is no general requirement of this. Also, there is 
no need to restrict to a single-component inflaton. From 
the particle theory viewpoint, all scalar fields are com- 
plex in supersymmetry, and there are many scalar fields. 
From the inflationary model building point of view, it 



C. Possible undesirable consequences of the 
traditional approach 

If \n'^ ln(A;/fc*)| > In* — 1| then the truncated Taylor se- 
ries gives a very unnatural form for the spectrum which 
is not motivated by any model of inflation, and which 
can give misleading parameter estimations. In general, 
the running can become appreciable for high k even if 
it is negligible for small k due to the possible k depen- 
dence of the running. Thus, if we ignore the running 
of the running, as is done in the traditional truncated 
Taylor expansion, the running of the spectrum at high 
k may be too biased by the data at low k because the 
running at high k is forced to be same as that at low 
k. In those cases, we should take account of the running 
of the running and more generally all the higher order 
terms. Considering infinitely many terms or an infinite 
number of parameters would, however, be unpractical 
in actual parameter estimations, and we shall suggest a 
more appropriate way of parameterizing the power spec- 
trum in Sec. IIIII Thus taking only the first few terms in 
the Taylor expansion could give a poor representation of 
the power spectrum and lead us to incorrect parameter 
estimations. 



In multi-component models, there are many slow-roll parameters. 
The relevant slow-roll parameters are the ones that directly affect 
the spectrum, see Ref. |3- 
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III. A BETTER APPROACH 

We suggest a different way to parameterize the spec- 
trum, which does not require any extra parameters com- 
pared with the traditional truncated Taylor expansion 
parameterization. In the standard slow-roll limit it re- 
duces to the usual truncated Taylor series form, but has 
a more sensible extension beyond the standard slow-roll 
limit. 



A. Improved parameterization 

Instead of the truncated Taylor expansion, we param- 
eterize the spectrum as 



lnV(k) = lnV* 



(n* - l) 2 



and hence the spectral index as 



n — 1 = (ra* — 1) 
and the running as 



k \ " 



(3.1) 



(3.2) 



(3.3) 



where, as before, fc» is an arbitrary reference point. Note 
that our parameterization has a simple form and uses 
the same number of parameters as the truncated Tay- 
lor series of Eq. (|2.1[1 . Expanding Eq. I|3.1|l in the limit 
|r?4 hi(k/k^)\ -C |n* — 1|, we see our parameterization re- 
produces the standard slow-roll Taylor series 

InV = \nP* + (n„ - 1) In + ~< In 2 (^y 

(3.4) 

of Eq. p.lfl . but with a more sensible extension to the 
domain |rt' t ln(fc/fc»)| > n* — 1|. 



B. Inflationary motivation of our parameterization 

The truncated Taylor series and our parameterization 
are equivalent for |n'„ ln(fc//c*)| <C \n* — 1|, but our pa- 
rameterization is also well motivated for |n'„ ln(A;/fc»)| > 
\n* — 1|- 

Specifically, in the general slow-roll approximation the 
spectrum for multi-component inflation models is given 
byi 2 



In V 



-ki W(k£)] 



hi IT 



2n' 

3 IF 



(3.5) 



where II = II(ln£), II' = dil/dln^, and £ is minus the 
conformal time: £ = — J dt/a ~ 1/aH where a is the 
scale factor and H is the Hubble parameter. The window 
function —xW'(x) is given in Refs. 0, @ and has the 
properties 



and 



r°° dx 

/ —\- x W\x)\ = l 
Jo x 



lim [-xW'{x)\ = 0(x 2 ) 



(3.6) 



(3.7) 



II represents the relevant inflationary parameters that 
directly affect the spectrum and is defined in Ref. 0. 

At zeroth order, II 2 can be regarded as constant and 
the normalization property of the window function leads 
to 



\nV = lnff 



(3.8) 



Our parameterization of Eq. (|3.1|) arises from a II 2 of 
the form 



inn 2 = inn; 



(3.9) 



where n 2 ^ is a constant and the second term with the 
constant coefficient B is assumed to be small, or 



n 2 s 

For v < 2, substituting into Eq. I|3.5|l gives 
ln7' = lnn 2 -Ak v 



(3.10) 



(3.11) 



where A is a constant depending linearly on B and non- 
trivially on v. For v > 2, the late time part of the integral 
dominates and n — 1 becomes proportional to k 2 . See 
Ref. [4| for a more detailed discussion. Comparing with 
Eq. (13. If) , we have 



— 1 



v for v < 2 
2 for v > 2 



(3.12) 



with standard slow-roll corresponding to v <C 1. Note 
that simple single-component inflation models tend to 



satisfy |n'| 



1| *C 1 7i — 1|. A concrete example of a 



particle theory motivated inflationary model which gives 
a spectrum of the form of Eq. I|3.12|l is given in Ref. |4| . 

Of course, the general slow-roll approximation can also 
accommodate cases where n'/n cannot be expressed as 
a power of £. In other words, this power law case is still a 
special case of the more general class of inflation models 
which the general slow-roll approximation can handle. 



IV. LIKELIHOOD ANALYSIS 



2 Note that in some cases there may be extra terms. See Ref. 
for the full story. 



We perform an estimation of cosmological parameters 
using different parameterizations of the primordial per- 
turbation spectrum: (1) our spectral parameterization 
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0.9 0.95 1 1.05 

FIG. 1: Shown are the likelihood contours of 68.3% and 95.4% 
for (n» ,n' t ) in our parameterization, using observations of the 
CMB, the SDSS 3D power spectrum of galaxies, and the mat- 
ter power spectrum inferred from the SDSS Lyman-a for- 
est. Also shown are the lines of n*/(n» — 1) = 0.2,0.5,1,2 
as solid, dashed, dot-dashed, and dotted, respectively. The 
curves n'»/(n» — l) 2 = ±1 are triple-dot dashed. 
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FIG. 2: Plotted are the marginalized one-dimensional proba- 
bility distribution functions for cosmological parameter esti- 
mation using our parameterization (solid), the truncated Tay- 
lor expansion (dash), and constant spectral index (dotted). 




FIG. 3: In the left panel, we show the marginalized one- 
dimensional probability distribution functions for n* for our 
parameterization (solid), for the truncated Taylor expansion 
(dash), and for constant spectral index (dotted). The right 
panel is that for n*/(n» — 1). 

given by Eq. JH2J; (2) the truncated Taylor expansion 
given by Eq. J23J; and (3) the case of constant spec- 
tral index. We model a flat universe with the cosmo- 
logical parameters fifc/i 2 , fl c h 2 , O s , hi^4 and r. Here, 
A is related to the amplitude of curvature perturbations 
at horizon crossing, |A H | 2 = 2.95 x 10~ 9 A at the scale 
= 0.05/i/Mpc. The angular acoustic peak scale 9 S is 
the ratio of the sound horizon at last scattering to that of 
the angular diameter distance to the surface of last scat- 
tering, and is a useful proxy for the Hubble parameter 
H = lOOh km s" 1 Mpc 0. 

We use a Markov- Chain Monte-Carlo (MCMC) tech- 
nique with a modified form of CosmoMC [8( • We use the 
CMB data from first year WMAP 9], ACBAR 0, CBI 
Jll| and VSA 0, the galaxy power spectrum from SDSS 
|13| and Lyman— a 14]. Each MCMC analysis included 
approximately 10 5 points, with an acceptance efficiency 
of approximately 45%. The parameter space was sam- 
pled from flat priors 
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which are all well outside of the regions of high prob- 
ability, except for the requirement of r < 0.3 for our 
parameterization case (1), which we include to exclude 
unphysically high optical depths otherwise allowable by 
the data and this form of the power spectrum. We use 
t < 0.8, outside appreciable levels of the PDF for cases 
(2) and (3). Since currently there are only upper limits 
on the contribution of tensor perturbations, we do not 
include tensor perturbations in our analysis. 

The 2-D likelihood of the parameters n'„ and rt* is 
shown in Fig. 1, where we also show lines of constant 
n'^/(n* — 1) and the curves of n^/(n^, — l) 2 = ±1. Note 
that in simple single-component inflation models one of- 
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TABLE I: Central values, 68.3% uncertainties, and \ 2 f° r 
the different models of the primordial scalar spectrum. 
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ten gets |n 



1| 2 <C |n — 1|, and our parameteriza- 



tion indicates the current data is consistent with a range 
well beyond that where this relation holds. In Fig. 2, we 
show the cosmological parameter estimation using the 
three different parameterizations, and in Fig. 3, we show 
the estimation of n* and n^/(n* — 1). The central values 
for the models' parameters, their uncertainty and % 2 are 
given in Tabled 

In Fig. 1, the 95.4% CL contour region has a stretched 
region towards positive n* — 1 and negative n'^ compared 
with the elliptic shaped 68.3% CL contour. This arises 
due to the combination of the fit trying to satisfy the low- 
l multipolcs simultaneously with the high-optical depth 
allowed. 

The data used in our analysis covers a k range of 
Aln/c ~ 10. This means the central values in Table I 
giving |n*| ~ 0.2 |n* — 1| for our parameterization, or 
\n'J ~ 0.4 |n* — 1| for the truncated Taylor series, are in- 
consistent with the range of validity of the truncated Tay- 
lor series, as they give |n'„ ln(fc/fc*)| ~ \n„ — 1|. Note that 
the central values of differ because, for n'^j (n*— 1) > 0, 
our parameterization alters the power spectrum more 
strongly for k > than the truncated Taylor expan- 
sion. Despite the lack of robust data at k even larger 
than the SDSS Lyman-a data, this result is a good il- 
lustration of the danger of ignoring the higher derivative 
terms in the Taylor series when we deal with more precise 
cosmological data covering a wide range of scales in the 
future. 



V. CONCLUSION AND DISCUSSION 

We discussed the possible bias and inconsistency in the 
cosmological parameter estimation induced by presuming 
a truncated Taylor series form for the power spectrum. 
We proposed an improved form for the power spectrum 



which is motivated from actual perturbation calculations 
applicable to a wide class of well motivated inflation mod- 
els. The standard slow-roll truncated Taylor series form 
is just a special limit of our parameterization of the power 
spectrum. Our proposed form requires no additional free 
parameters compared with the traditional truncated Tay- 
lor expansion, and can be straightforwardly implemented 
in existing codes. Our form for the spectrum is simple 



In 7? = ln"P - Ak v 
= -vAk v 



(5.1) 
(5.2) 



with only three constant free parameters, Vq, A and 
v. There exist several possible ways to parameterize 
this form. The parameterization used in this paper, 
Eq. (|3.1|l . was motivated by clarity in the comparison 
with the truncated Taylor expansion. The standard slow- 
roll cases correspond to v — n^/(n* - 1) « 1. In per- 
forming the likelihood analysis via MCMC, we found 
the central values for n^/(n* — 1) were inconsistent with 
|n.* ln(fc/fc*)| <S |n* — 1|, the basic assumption behind the 
truncated Taylor series. 

Further to the technical discussion on the form of the 
power spectrum using the general slow-roll formula in 
Sec. IIII Bl let us point out that the form we presented 
in this paper should be regarded as an asymptotic form 
for small Ak v . This asymptotic form is sufficient for the 
range of our analysis in this paper, i.e. the data up to 
k ~ 6h/Mpc from Lyman— a. If we can include data 
at higher k, well beyond the Lyman-a range, covering a 
change in In 7-* of order unity, we may need to use the 
original form of Ref. Q without asymptotic approxima- 
tion. Otherwise V would decrease too rapidly for higher 
k (assuming v > Q). Alternatively we may soften the 
form of our parameterization, for example as 



lnT 3 = ln"P Q 

Avk v 

n — 1 = — - 



(5.3) 
(5.4) 



1 + Ak v 

so that it gives Eqs. (|5.1() and (|5.2() for small Ak v ' , and 



In 7? = lnT 3 * 
n — 1 = — v 



v In 



(5.5) 
(5.6) 



for large Ah", or ideally add more parameters. 

Simple single-component inflation models require the 
standard slow-roll approximation to produce a flat spec- 
trum, but this is not the case for multi-component infla- 
tion models. In this sense, our parameterization would 
be of great interest for multi-component inflation models, 
such as given in Ref. p|. 

Our parameterization, the traditional truncated Taylor 
expansion, and the case of constant spectral index led to 
quite similar, though not identical, results for the cosmo- 
logical parameter estimations using the currently avail- 
able data dominated by large scale observations such as 
CMB and galaxy surveys. This indicates that currently 
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the running is not crucial to the cosmological parameter 
estimations, but this may change with reduced error bars 
in the near future 0] . Application of a parameterization 
such as ours would also be important for modeling very 
small scale structure, such as the first objects in universe 
[Tfil | which could extend the range up to k ~ 10 6 h/Mpc. 
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